Branching Ratio of the Electromagnetic Decay of the S+(1385) 



D. Keller, K. Hicks, K.P. Adhikari,^^ D. Adikaram,^* M.J. Amaryan,^^ M. Anghinolfi,^* H. Baglidasaryan,^®'28 
J. Ball,6 M. Battaglieri,i8 I. Bedlinskiy,^! A.S. Biselli,!^'^ c. Bookwalter,!^ S. Boiarinov,^^ D. Branford,^ 

W.J. Briscoe," W.K. Brooks,36.34 y.D. Burkcrt,^* D.S. Carman,^^ A. Celentano,!^ S. Cliandavar,^^ 
P.L. Cole,i5'34 Contalbrigo,!'^ V. Crede,!^ A. D'Angelo,!'''^! a. Danicl,^^ N. Dasliyan,40 R. De Vita,!^ 
E. De Sanctis,!'^ C. Djalali,^^ D. Doughty,^- R. Dupre,i A. El Alaoui,i L. El Fassi,i L. Elouadrliiri,^-* 
P. Eugenio,i2 G. Fedotov,33 M.Y. Gabrielyan," N. Gevorgyan,40 G.P. Gilfoyle,^^ K.L. Giovanetti,^^ W. Gohn,* 
E. Golovatcli,32 R.W. Gotlie,33 L. Graham,^^ K.A. Griffioen,^^ M. Guidal,^^ N. Guler,28.EI L. Guo,"'34 
K. Hafidi,! H. Hakobyan,36.40 ^ Holtrop,^^ Y. Ilieva,^^. i4 q ireland,^^ B.S. Isliklianov,32 E.L. Isupov,32 
H.S. Jo,20 K. Joo,® M. Khandaker,26 p. Khetarpal," A. Kim,23 W. Kim,23 F.J. Klein,^ A. Kubarovsky,29'32 

V. Kubarovsky,34.29 g y. Kuleshov,36. 21 h.Y. Lu,^ I .J .D. MacGregor 37 y. Mao,33 N. Markov,** 
M. Mayer, 28 B. McKinnon,^^ G.A. Meyer,'* M. Mirazita,^'^ V. Mokeev,34. 32,[t| jj^ Moutarde,*^ E. Munevar," 
p. Nadel-Turonski,34 R. Nasseripour,"' " S. Niccolai,20 G. Niculescu,22^ 27 i Niculescu,22' 34, i4 Qsipenko,!® 
A.I. Ostrovidov,i2 y[_ Paolone^^ L. Pappalardo,*'^ R. Paremuzyan,'*° K. Park,34.23 Park,i2 Pasyuk,34.2 
Anefalos Pereira,^^ S. Pisano,2°'F| 0. Pogorelko,2i S. Pozdniakov,2i S. Procureur,^ Y. Prok,^'38 Protopopescu,'^'' 
B.A. Raue,"'34 q Ricco,"^irD. Rimal," M. Ripani,i^ B.G. Ritcliie,^ G. Rosner,^^ p. Rossi," F. Sabatie,^ 
M.S. Saini,i2 C. Salgado,^^ D. Schott," R.A. Schumacher,* H. Seraydaryan,28 Y.G. Sharabian,^* E.S. Smith,^* 

G.D. Smith,37 D.I. Sober,5 D. Sokhan,20 S.S. Stepanyan,23 S. Stepanyan,34 P. Stoler,29 S. Strauch,33' 14 
M. Taiuti,i3.| W. Tang,27 C.E. Taylor,!^ S. Tkachenko,^^ B .Vernarsky,* M.F. Vineyard,^^. 30 a.V. Vlassov,2i 
H. Voskanyan,*".!!! E. Voutier,24 M.H. Wood,3.33 n. Zachariou,** L. Zana,25 B. Zhao,39 and Z.W. Zhao3S 

(The CLAS Collaboration) 

^ Argonne National Laboratory, Argonne, Illinois 60439 
^Arizona State University, Tempe, Arizona 85287-1504 
^ Canisius College, Buffalo, NY 
Carnegie Mellon University, Pittsburgh, Pennsylvania 15213 
^Catholic University of America, Washington, D.C. 20064 
"CEA, Centre de Saclay, Irfu/Service de Physique Nucleaire, 91191 Gif-sur-Yvette, France 
^Christopher Newport University, Newport News, Virginia 23606 
^University of Connecticut, Storrs, Connecticut 06269 
'^Edinburgh University, Edinburgh EH9 3JZ, United Kingdom 
"'Fairfield University, Fairfield CT 06824 
Florida International University, Miami, Florida 33199 
^^Florida State University, Tallahassee, Florida 32306 
^'^ Universita di Cenova, I6I46 Cenova, Italy 
'-^ The Ceorge Washington University, Washington, DC 20052 
^^Idaho State University, Pocatello, Idaho 83209 
^^INFN, Sezione di Ferrara, 44100 Ferrara, Italy 
^"^INFN, Laboratori Nazionali di Frascati, 00044 Frascati, Italy 
^^INFN, Sezione di Cenova, I6I46 Cenova, Italy 
'^INFN, Sezione di Roma Tor Vergata, 00133 Rome, Italy 
'"^"institut de Physique Nucleaire ORSAY, Orsay, France 
^''institute of Theoretical and Experimental Physics, Moscow, 117259, Russia 
^'"^ James Madison University, Harrisonburg, Virginia 22807 
^'^ Kyungpook National University, Daegu 702-701, Republic of Korea 
''^^LPSC, Universite Joseph Fourier, CNRS/IN2P3, INPC, Grenoble, France 
University of New Hampshire, Durham, New Hampshire 03824-3568 
'''^Norfolk State University, Norfolk, Virginia 23504 
^^Ohio University, Athens, Ohio 4 5701 
Old Dominion University, Norfolk, Virginia 23529 
Rensselaer Polytechnic Institute, Troy, New York 12180-3590 
""^University of Richmond, Richmond, Virginia 23173 
^' Universita' di Roma Tor Vergata, 00133 Rome Italy 
Skobeltsyn Nuclear Physics Institute, Skobeltsyn Nuclear Physics Institute, 119899 Moscow, Russia 
''^University of South Carolina, Columbia, South Carolina 29208 
Thomas Jefferson National Accelerator Facility, Newport News, Virginia 23606 

Union College, Schenectady, NY 12308 
""^ Universidad Tecnica Federico Santa Maria, Casilla 110-V Valparaiso, Chile 



2 



University of Glasgow, Glasgow G12 8QQ, United Kingdom 
''^^ University of Virginia, Charlottesville, Virginia 22901 
College of William and Mary, Williamsburg, Virginia 23187-8795 
'^''Yerevan Physics Institute, 375036 Yerevan, Armenia 
(Dated: November 28, 2011) 

The CLAS detector was used to obtain the first ever measurement of the electromagnetic decay 
of the E*+(1385) from the reaction 7p — >■ isr''E*+(1385). A real photon beam with a maximum 
energy of 3.8 GeV was incident on a liquid-hydrogen target, resulting in the photoproduction of 
the kaon and E* hyperon. Kinematic fitting was used to separate the reaction channel from the 
background processes. The fitting algorithm exploited a new method to kinematically fit neutrons 
in the CLAS detector, leading to the partial width measurement of 250.0 ± 56.9(staf)t4i'2 (sys) 
keV. A U-spin symmetry test using the SU(3) flavor-multiplet representation yields predictions 
for the E*+(1385) E+7 and E*''(1385) A7 partial widths that agree with the experimental 
measurements. 

PACS numbers: 13.40.Em,14.20. Jn,13.30.Ce,13.40.Hq 



INTRODUCTION 

The electromagnetic (EM) decay of baryons can pro- 
vide considerable information on their underlying struc- 
ture. This transition offers a clean probe of the wave- 
function of the initial and final state baryons, provid- 
ing theoretical constraints and tests of the quark model. 
The non-relativistic quark model (NRQM) of Isgur and 
Karl ffl, f^l predicts the electromagnetic properties of the 
ground state baryons reasonably well. It has been less 
successful giving accurate descriptions of the low-lying 
excited-state hyperons. Several other theoretical tech- 
niques have been used to more accurately calculate these 
transitions, including (NRQM) jES], a relativized con- 
stituent quark model (RCQM) j5!|, a chiral constituent 
quark model (xCQM) |6|], the MIT bag model 7], the 
bound-state soliton model a three- flavor generaliza- 
tion of the Skyrme model that uses the collective ap- 
proach!^, [1^ , and an algebraic model of hadron struc- 
ture [nr 

Photoproduction from nucleon targets is a useful tech- 
nique to cleanly generate a significant statistical sample 
of hyperons and to measure EM transitions to other de- 
cuplet baryons. If the EM transition form factors for 
decuplet baryons with strangeness are also sensitive to 
meson cloud effects, models attempting to make predic- 
tions of the decuplet radiative decay widths will need re- 
visions to incorporate this effect. Comparison of data for 
the EM decay of decuplet hyperons, S*, to the present 
predictions of quark models provides a measure of the 
importance of meson cloud diagrams in the E* — 
transition. Experimental results for the EM decay ratios 
for all E* charge states are desirable to obtain a complete 
comparison to EM decay predictions for the E*. Preci- 
sion measurements of the E*^ — >■ E^7 and E*+ E+7 
decay widths can be particularly useful in determining 
the degree of SU(3) symmetry breaking. 

The decay width from the measurement of E**^ — > A7 
l3, l3| is much larger than most current theoretical pre- 
dictions. This could be due to meson cloud effects, which 



were not included in these calculations. There is a the- 
oretical basis for calculating these effects [3] that sug- 
gests pion cloud effects may be sizable. For example, 
they are predicted to contribute on the order of '^40% 
to the 7p — ^ A+ magnetic dipole transition form fac- 
tor, G^jIq^), for low Q2, The CQM [T^ indicates that 
the value of Gfj{Q) is directly proportional to the proton 
magnetic moment , and measurements of G^j for low 
are rationalized in the framework of the model if the 
experimental magnetic moment is lowered by about 25%. 

With theoretical predictions for the degree at which 
the meson cloud effects plays a role, it is then possible 
to test SU{3) flavor symmetry breaking (and the degree 
at which it is broken). This can be achieved by measur- 
ing both the r(E*- ->■ E-7) and r(E*+ E+7) decay 
widths and comparing these to predictions from flavor 
SU(3) relations. 

Just like isospin invariance can be used to compare the 
A++ p7r+ and A+ — ^ p7r° decays, U-spin invariance 
may be used to compare the E*+ — > E+7 and A+ — > pj 
decays. U-spin is analogous to isospin in that it is a sym- 
metry in the exchange of the d and s quarks rather than 
the u and d quarks. A value of U-spin can be assigned 
to each baryon based on its quark composition. The 
E*^ and the S*~ of the baryon decuplet have U = 3/2, 
whereas the octet baryons E^ and have U = 1/2. 

U-spin symmetry forbids radiative decays of specific 
decuplet baryons. Since the photon is a charge singlet 
with U = 0, this implies that 



E*" ^ E"7 and S* 
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have zero amplitude in the equal-mass limit due to U-spin 
symmetry. This can also be understood in the context 
of the SU(6) wavefunctions for these baryons. The Ml 
transition operator is written between the initial and final 
states as : 
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Here the sum is over all q constituent quarks, rUq, Oq and 
Qq are the mass, spin vector and charge of the q quark, 
k is the propagation direction, and e*"^ is the polariza- 
tion vector. One can also show that the same transition 
operator for the S*"*" gives a non-zero amplitude. U-spin 
invariance implies a large difference in the radiative decay 
widths of the E*" and !]*+. 

The chiral symmetry for U-spin is strongly broken be- 
cause the constituent mass of the strange quark, rris, is 
approximately 1.5 times greater than the non-strange 
quarks, m. The magnetic moment is inversely propor- 
tional to the mass, and so there is no cancellation in the 
wavefunction like in the equal-mass SU (6) case in Eq. 
[TJ From Ref. [31, an estimate of the ratio of the EM 
decay rates from the ratio of the square of the transition 
operators can be expressed as 
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resulting in a value of about 1%. This suggests that U- 
spin symmetry breaking for radiative decays is at the 
level of only a few percent. At this level U-spin is an 
effective tool, even considering the quark mass difference. 

Detailed calculations from the CQM and l/7Vc-type ex- 
pansions of the EM decay rates have been carried out by 
several groups 17|, |l8|, all of which come up with de- 
cay ratios of a similar scale. In lattice QCD, the quarks 
have very different interactions with the photon than for 
the CQM, but these too have ratios (for the above equa- 
tion) within a few percent [l^. This consistency makes 
a stronger case for the usefulness of U-spin symmetry. 

There has been much theoretical interest in radiative 
baryon decays. However, there are only a few measure- 
ments. Recently, a measurement of the radiative decay 
of the S*~ was attempted by the SELEX collaboration 
[20| . resulting in only an upper limit. The 90% confi- 
dence level upper bound of F = 9.5 keV was reported, 
however most models predict a value of less than 4 keV. 
Ultimately this result has limited power to constrain the- 
oretical estimates. More experimental measurements are 
necessary to provide better constraints. 

A program to investigate the various S* electromag- 
netic decays is underway using data from the CEBAF 
Large Angle Spectrometer (CLAS) detector. First, two 
independent analyses of the EM decay of the S*" have 
been completed 12 1, [13] . The consistency in these re- 
sults has given confidence in the notion that meson cloud 
effects are indeed contributing significantly. The next 
step described and presented here was to measure the 
S*+ electromagnetic decay, which has not been done be- 
fore. The final program analysis for the E*~ — >■ E~7 
decay will be addressed in a future CLAS publication. 

In the following, a description of the experimental de- 
tails and analysis procedure for extracting the E*"*" EM 
decay branching ratio normalized to the strong decay is 
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FIG. 1: The CLAS detector at Jefferson Lab showing 
the toroidal magnet, the drift chambers, the time-of-flight 
scintillators, the cerenkov counter, and the electromagnetic 
calorimeter. 



provided. Some specifics are given about neutron detec- 
tion and the development of the neutron covariance ma- 
trix required by the analysis. After the signal extraction 
a U-spin symmetry test using the U-spin SU(3) multi- 
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plet representation is used to predict the E 
and E*'' A7 partial widths, which are then compared 
to the experimental results. 



THE EXPERIMENT 

The present measurements were carried out with the 
CLAS in Hall B at the Thomas Jefferson National Ac- 
celerator Facility [2l|. An electron beam of energy 4.023 
GeV was used to produce a photon beam with an energy 
range of 1.6-3.8 GeV, as deduced by a magnetic spec- 
trometer [2§| that "tagged" the electron with an energy 
resolution of ^ 0.1%. A 40-cm-long liquid hydrogen tar- 
get was placed such that the center of the target was 10 
cm upstream from the center of CLAS. 

The CLAS detector is constructed around six super- 
conducting coils that generate a toroidal magnetic field 
to momentum-analyze charged particles. The detection 
system consists of multiple layers of drift chambers to 
determine charged-particle trajectories, Cerenkov detec- 
tors for electron/pion separation, scintillation counters 
for flight-time measurements, and calorimeters to iden- 
tify electrons and high-energy neutral particles, see Fig. 
[H The Cerenkov detectors are not required for this ex- 
periment. 

Each event trigger required a coincidence between the 
OR of the detector elements in the focal plane of the 
photon spectrometer and the CLAS Level 1 trigger. The 
Level 1 trigger required two charged particles in two dif- 
ferent sectors of CLAS within a 150 ns coincidence time 
window. The approximate integrated luminosity for the 
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CLAS glla run period used in this analysis was 70 pb ^. 
Details of the experimental setup can be found in Refs. 



21, 2 



EVENT SELECTION 

Events were selected for the channel — > 
The present Particle Data Group (PDG) branching ra- 
tios list the decay S* Stt to be 11.7 ± 1.5%, and as- 
suming isospin symmetry, this leads to a branching ratio 
of 5.85 ± 0.75% for the £*+ ^ T.+Tr° decay This 
channel will be used to normalize the radiative signal that 
comes from the channel E*"*" S+7. For both channels 
the topology of the decay is jp — > K'^Y^^ (X), where X 
is not directly measured, such that the tt'^ and 7 are dif- 
ferentiated using conservation of energy and momentum. 
This topology leads to the final set of decay products 
7P K°i:+lx) 7r+7r-7r+n(X). The charged parti- 
cles can easily be detected with the use of the CLAS drift 
chambers and timc-of-flight system. The neutron must 
be detected with the CLAS electromagnetic calorimeters. 
The analysis was done using a previously prepared data 
reduction (skim) that required two positively charged 
tracks and one negatively charged track for each event. 

Cuts were applied to take into account both the re- 
gions of CLAS where there are holes in the acceptance 
that arose from problematic detector elements and re- 
gions that were not well simulated. This includes tracks 
at extremely forward or backward angles, areas near the 
torus coils, and regions where the drift chambers and 
scintillator counter efficiencies were not well understood. 
Tracks that point near these shadow regions are less likely 
to be reconstructed accurately. In addition a minimum 
momentum of 0.125 GeV, after energy loss corrections, 
was enforced for both positively and negatively charged 
particles to ensure accurate drift chamber track recon- 
struction. 

During the initial data skim, the hit times in the start 
counter that surround the target were used to find an 
interaction vertex time for each charged particle, which 
was then matched up with photons identified in the tag- 
ger, where there can be up to 10 candidate photons for 
a given event. The photon with the closest time to any 
track was selected as the photon that caused the event. 
Specifically, the time of interaction was determined using 
the time of the electron beam bucket (the accelerator RF 
time) that produced the event. To correlate the inter- 
action time with the photon production time, a timing 
coincidence between the tagger and the start counter was 
used. The RF time for the photon was then used to get 
the vertex time (photon interaction time tj) for the event. 
Using the time-of-flight (TOF) from the event vertex to 
the scintillator counter, the velocity /3 was calculated for 
each particle. From /3 and the particle's measured mo- 
mentum, a mass was calculated. Each track did not need 



to have a hit registered in the start counter for its mass 
to be calculated, only one track in the event needed a 
start counter hit. 

The mass squared calculated from time-of-fiight is 



''cal 



(2) 



where /? = L/tmeas such that L is the path length from 
the target to the scintillator, t„ieas is the measured time- 
of-flight, and the speed of light is set to 1. From this ini- 
tial identification, it was possible to use additional timing 
information to improve event selection. The measured 
time-of-flight and calculated timc-of-flight were used for 
an additional constraint. The measured time-of-flight is 
trneas = tgc — tj, whcrc tsc is the time at which the par- 
ticle strikes the time-of-flight scintillator counter. At is 
then 



At = tr. 



t. 
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where tcai is the time-of-flight calculated for an assumed 
mass such that 



teal = L\ 1 



(4) 



where m is the assumed mass for the particle of interest, 
and p is the momentum magnitude. Cutting on At or 
rricai should be effectively equivalent. 

Using At for each particle it was possible to reject 
events that were not associated with the correct RF beam 
bucket, which was separated by 2 ns. This was done by 
requiring \At\ <1 ns for all charged particles in the ini- 
tial analysis. This cut was chosen to minimize signal loss 
while also minimizing overlap from other beam buckets. 

A A/3 cut was used to clean up the identification 
scheme. A/3 is the difference between the measured 
/3 — L/{tmeas) and the calculated P — p/ \/p^ + rn? ■ The 
good events were taken within a cut of —0.035 < A/3 < 
0.035 for all pions as shown in Fig. [21 



Kaon identification 

In the reaction of interest, 7p — >■ i(r°S*+(X) 
7r+7r^7r+n(X), it is necessary to determine which tt"*' 
comes from the . It is possible to check both final 
state 7r"'''s with the detected -k" to study the kaon can- 
didates in each case by using the invariant mass. 

The invariant mass was selected for each tt+tt^ pair, 
as shown in Fig. [31 Whichever 7r+ lead to the in- 
variant mass that was closest to the mass of the 
was associated with the identification. Afterwards 
a cut at ±0.01 GeV about the mass was used to 
clean up the selection. For cases where both tt'^ com- 
binations with the 7r~ fell within the K'^ mass limit the 
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FIG. 2: The A/3 distributions for 7r+ (left) and tt" (right). 
The cut of —0.035 < A/3 < 0.035 is shown as the dashed lined 
in each case. 

wrong 7r+ could be selected. Monte Carlo was used to 
check the frequency of this ambiguity and was found 
to be ~ 2% for the 7p K°J:*+ 7r+7r-7r+n7r° and 
7P — >■ if°E*+ 7r+7r~7r+n7 channels. With additional 
kinematic constraints these ambiguous events were ulti- 
mately rejected. 




M(T:*r) (GeV) 

FIG. 3: Invariant mass of the ■k'^--k~ combination for the two 
different tt^ detected, prior to any tt"*" organization. 



Neutron identification 

Neutral particles are detected in CLAS as clusters 
in the electromagnetic calorimeter (EC) [lH not associ- 
ated with any reconstructed charged track from the drift 
chambers. The momentum reconstruction depends on 
the path length and TOF of the neutron. The direc- 



tional components of the neutral track were found by 
using the vertex and the cluster position on the EC for 
that hit. In this experiment the information about the 
neutral vertex was limited to the information that could 
be extracted from the other charged particle vertices in 
the decay chain. 

The EC has six triangular sectors made of alternating 
layers of lead and scintillator. Scintillator layers compose 
of about 10-cm-wide scintillator strips, where strips in 
every consecutive layer run parallel to one of the three 
sides of the triangle. The EC has 13 layers of scintillator 
strips for each of the three directions making 39 total 
layers. In each direction the EC is subdivided into an 
inner stack of 5 layers and an outer stack of 8 layers. 

The EC reconstruction software forms a cluster by first 
identifying a collection of strips in each of the three views. 
The software requires a set of threshold conditions to be 
met and that the strips to be contiguous. The groups of 
strips that pass these conditions define a peak and are 
organized with respect to the sum of the strip energies. 
The peak centroid and RMS in each of the three views 
is obtained and clusters are identified as intersection of 
centroids of peaks within their RMS. If a given peak con- 
tributes to multiple hits, then the energy in each hit due 
to that peak is calculated as being proportional to the rel- 
ative sizes of the multiple hits as measured in the other 
views. For example, if there are multiple hits which have 
the same U peak, the energy in V and W is added for 
each of the hits, and the ratio of these summed ener- 
gies determines the weight of the U peak's energy of the 
multiple hits. If the software thresholds for the scintilla- 
tor strip, peak and weighted hit energy are met then the 
cluster position and time are recorded. The events EC 
time (or EC time-of-flight) is defined as the time between 
the event start time and the time of the EC cluster. 

During analysis the strip information was used to de- 
termine whether the centroid was reconstructed using 
only the outer stack of the EC or both the inner stack 
and outer stack. The centroid could be located in any 
one of the layers of each stack, however, the cluster re- 
construction position did not contain that information, 
so the hit was assumed to be on the upstream face (closer 
to the target) of whichever stack the hit was contained 
in. With the assumed reaction vertex and the EC clus- 
ter position, the directional components in 9 and (j> were 
found, as well as the path length of the neutron. Us- 
ing the EC time-of-flight the momentum was calculated. 
The neutrons were differentiated from photons using a 
/3 < 0.9 cut. 

Neutron detection is essential for the reaction of in- 
terest. The neutron momentum was used in combina- 
tion with the 7r+ not associated with the to study 
the kinematics of the E+. Having clean constraints on 
the and is important when considering the event 
topology 7p ^ K°E*+ -)■ K°Y;+{X). 

A thorough study of the accuracy of the EC for neu- 
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tron reconstruction in all kinematic ranges has not been 
achieved previously at CLAS. Obtaining the resolution 
in all measured variables for neutron reconstruction was 
an essential part of the present analysis. Correlations 
between each measured variable in the EC had also not 
been previously studied. The EC covariance matrix of 
the neutron can give a lot of information about the qual- 
ity of the kinematic variables in each region of the EC. 
These values can then be used to weight the neutron mea- 
surements appropriately in kinematic constraints that de- 
pend on maximum likelihood methods [26]. 

There are resolution differentials in all measured vari- 
ables that are related to the acceptance of the EC. Hits 
from the center of each triangular sector have better mea- 
surements over those on the edges due to shower leakage. 
The inner and outer stacks can act as separate detec- 
tors in the sense that if a hit is seen in the outer but 
not the inner stack, then the inner stack plays no role 
in the reconstruction of that hit. It is far less common 
for an event to pass through the inner stack with no ef- 
fect and to register a hit in the outer stack, but for these 
events, the outer EC stack was used independently with 
its own unique resolution parameters for each measured 
variable. All possible combinations of the measured neu- 
tron dependence on 0, 0, and p were studied to develop 
a complete understanding of the neutron variance and 
covariance in the EC 12711. 



Neutron detection test 

The test reaction — ^ 7r+7r^7r+n was isolated in the 
glla data set by selecting a tt^ and two 7r+ , and kinemat- 
ically fitting to a missing neutron hypothesis and then 
taking a 10% confidence level cut. Only the detected 
neutrons found in a direction less than 3° from the kine- 
matically fit three pion missing momentum were used to 
ensure the correct neutron. This channel was selected 
because the final decay products are identical to the re- 
action of interest jp K°J:*+ iC°S]+(X). In ad- 
dition, the momentum range of the detected particles is 
the same. Kinematic constraints were imposed to remove 
possible TT"*" 7r~ combinations with invariant mass equal 
to the K'^ , so that only the jp — ^ 7r"'"7r~7r"'"n events sur- 
vived. The simplification made by working with the test 
channel is that in the "fp — n reaction, there is 
only one interaction vertex. This implies that the neu- 
tron comes from the primary interaction vertex, which 
can be well determined using the charged pions. 

To study the measured neutron variable residuals, we 
required each event to have one detected neutron and 
then compared the measured variable with the kinemat- 
ically fit missing variable in each case. Assuming a high 
quality missing neutron four vector, this procedure was 
used to find the change in resolution with respect to all 
measured variables over the EC face ^] . Only the events 



that registered an actual hit in the EC were used to study 
the resolution. No EC fiducial cuts were applied during 
the covariance investigation so that the entire EC face 
could be studied and compared to Monte Carlo. During 
analysis, only the minimal fiducial cuts were applied of 
8° < 6* < 40° on the neutron polar angle to maximize the 
statistics. 

For the test channel the neutron vertex was found from 
a multi-track- vertex fitting procedure to give an accurate 
vertex (at less than 4% uncertainty in position for the 
topology of interest) for multi ple final state particles all 
coming from the same vertex [28| |. Because the neutron 
came from the primary interaction vertex in this study, 
its vertex was accurately known. However, for events in 
which the neutron comes from secondary vertices, its ver- 
tex is not as easily obtained. Because the neutron vertex 
information can affect its reconstructed four momentum, 
these differences can be important when studying reso- 
lutions. 

Once the EC neutron covariance matrix for jp — > 
TT+TT^TT+n was well understood, the Monte Carlo res- 
olution was matched to the data using the same test 



channel [27|. The Monte Carlo was then used to study 
the jp K°J:*+ K°^+TT° channel and to find the 
neutron covariance matrix specific for this topology. In 
this way the K'^ interaction point with the beam line 
could be used as the starting point of the neutron path 
in the neutron reconstruction process for any of the 
decays, so no bias was introduced by assuming a S"*". 
This step removed the explicit dependence on the neu- 
tron vertex. The Monte Carlo covariance matrix for 
7P K°Y.*+ iir"S+7r° was then used to tune the 
data neutron covariance matrix specific to the topology. 
The change in the momentum resolution from this tuning 
process was smaller than 5%. 

In order to obtain a consistent covariance matrix for 
the neutron, discrimination was made for each neutron 
between the inner and outer EC stacks in order to cal- 
culate the correct path length. In addition, timing and 
momentum corrections were applied as described below. 



Neutron path 

As previously stated, the distance that the neutron 
travels in CLAS was used with the EC time-of-flight to 
determine the momentum of the neutron. The distance is 
dependent on the EC stack and the position of the cluster 
reconstruction. The inner stack cluster reconstruction 
was always used unless there was only a hit in the outer 
stack. A determination of whether there was a hit only 
in the outer stack, the inner stack, or both, was made by 
checking which EC scintillators were associated with an 
event. The probability to find a hit in the outer stack 
alone was less than 15%. For all other neutral hits either 
the inner stack or both were associated with the hit. If 



7 



there was a hit only in the outer EC stack, the first layer 
(layer closest to the target) of the outer stack gave the 
plane of the EC cluster coordinates. If there was a hit in 
the inner stack or both, the first layer of the inner stack 
was used as the plane of the EC cluster coordinates. 

Neutron time 

The time-of-fiight for the neutron came from the dif- 
ference between the event start time and the EC cluster 
time. The path length used to reconstruct the neutron 
momentum, which assumes a hit on the EC face of either 
the inner or outer stack, was inaccurate by the distance 
the neutron traveled past the EC face into the detector. 
A correction was used to compensate for the average ad- 
ditional distance the neutron travels into the EC. In addi- 
tion the outer stack is farther from the target and for the 
same event would have a slightly different time response 
than that of the inner stack. 

A correction was implemented directly in the neu- 
tron time-of-flight to correct the neutron momentum. 
This was done by using the calculated neutron time-of- 
flight and comparing it to the expected time Texpected — 
L/cjSmiss- Here, L is the path length of the neutron, 
and Pmiss is the /3 calculated using the missing mo- 
mentum and energy of the neutron from the 773 — > 
7r"'"7r~7r+(n) events that passed a 10% confidence level 
cut from a kinematic fit under a missing neutron hy- 
pothesis. By using Texpected — Tec for each stack, a 
separate correction was found for each case, such that 
Texpected - {Tec + Tcorrection) ~ 0. By finding Separate 
timing corrections for the inner and outer stacks, the far- 
ther distance of the outer stack was compensated for. 
The separate study of timing corrections for the inner 
and outer stacks was carried out using the Monte Carlo 
simulations. 

The timing correction used is the same in all direc- 
tions. However, in order to obtain accurate covariance 
information an additional momentum correction was re- 
quired which is sensitive to the geometry of the EC and 
neutrons trajectory. It was only after all corrections that 
the residual means of all measured variables were cen- 
tered around zero to accurately refiect the neutron reso- 
lutions. 

Neutron momentum correction 

A neutron momentum correction was implemented by 
studying the trend found in momentum and position res- 
olutions over various kinematic ranges. This was done by 
studying the residuals Ap, A0, and A0 over each variable 
p, 9, and (f). The residual of momentum, Ap, is defined 
as the difference between the missing neutron momentum 
and the reconstructed neutron momentum. Likewise for 
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FIG. 4: Top left: the neutron Ap distribution before cor- 
rection, top right: the Gaussian mean from the binned fits, 
bottom left: the corrected distribution, and bottom right: the 
final corrected Gaussian mean from the binned fits. 



the directional components A9 and A(f). The missing 
neutron four-vector was found by kinematically fitting 
the charged decay products in the missing mass of the 
neutron and taking a 10% confidence level cut. In this 
fit there were three unknowns from the components of 
the missing momentum vector and four constraints from 
the conservation of energy and momentum to make a 
1-C fit .26]. The trend of each of the residuals should 
be distributed around zero, if it is not, the distribution 
will display a trend that can be used to correct the mea- 
sured variable. Once the neutron momentum magnitude 
and directional resolutions are evenly distributed around 
zero, the missing and detected four- vectors are compa- 
rable. This implies that for the majority of events, the 
detected neutron momentum vector was the same within 
the experimental resolution as the high quality kinemat- 
ically fit missing neutron momentum vector. 

Similar corrections were implemented to Ap with re- 
spect to (j> such that p' = p + /(</)), A0 with respect to 
(j> such that (p' = (f) + f{4>), and 9 with respect to 9 such 
that 9' ~ 9 + f{9). The Monte Carlo required separate 
corrections in the same variables that were determined 
using the same procedures as for the data. 
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Neutron covariance matrix 

The neutron covariance matrix was determined after 
the corrected neutron path was used with the timing cor- 
rection implemented for the corresponding EC stack and 
momentum correction. This covariance matrix was re- 
quired in order to kinematically fit the neutron with the 
other detected particles. The variables used to represent 
the neutron vector components were 0, (j), and p, leading 
to a covariance matrix of 

vr = vf vf' vf' . 

The variance and correlations of each variable were ob- 
tained by studying the differences between the kinematic 
variables of the detected neutron from the kinematically 
fit missing neutron. The residuals in each case were sliced 
and binned to fit with a Gaussian to find the functional 
dependence in each variable. Once the functional depen- 
dence on each variable was found for each (j) and 9 bin, an 
empirical smearing technique was used to get the Monte 
Carlo to closely match the same functional dependence 
seen in the data. Similar steps were taken for the direc- 
tional components (27i |. 

ANALYSIS PROCEDURE 

In the analysis, progressive steps were taken to remove 
as much identifiable background as possible while pre- 
serving the counts from the channels 7p iir°I]*+ 
is:°S+(7r°) and jp K°i:*+ ^ K°Y;+{-f). The radia- 
tive signal was buried by the tt" — 77 decay and required 
advanced fitting techniques to resolve the signal. The 
fitting procedure developed here required that all other 
backgrounds be removed or extensively minimized to en- 
sure high quality separation between the radiative and 
strong decays of the E*+. 

The goal was then to achieve clean hadron identifica- 
tion before using the fitting procedure for the competing 
7r° and radiative signals. For the sake of notation, let tt^ 
indicate the tt'^ used in the invariant mass selection, 
such that TT^ is the 7r+ that forms the closest known K'^ 
mass when combined with the tt^. Naturally tt^ is the 
other detected 7r+. Fig. [5] shows the invariant mass of 
the TTi-TT^ (upper left), missing mass off the n^-n^ (up- 
per right), the n-7r^ invariant mass (lower left), and the 
missing mass squared of all the detected particles (lower 
right) . The distributions in Fig. [S] are before any kine- 
matic constraints and after the tt"*" assignments are made. 
The K''^ was cut about ±0.01 GeV of the known K'^ mass 
to reduce the tt^-tt^ background. Fig. [6] shows the in- 
variant mass of the n-7r^ (dashed lines show the cut that 
was implemented) (upper left), missing mass off the 



(upper right), the missing energy off all detected particles 
(dashed lines show the cut that was implemented) (lower 
left), and the missing mass squared of all the detected 
particles (lower right), after the K° cut. The E"*" peak 
is clearly visible as seen in the (upper left) plot. The 
clear visible peak in the missing mass squared at the 7r° 
mass is also an indication that the neutron measurement 
is effective. 

A Monte Carlo study on the phase space of the jp ^■ 
^Oj^*+ _^ K'^'E+(^X) reaction indicated that most events 
from the missing energy boosted in the S*"*" frame 
should be in the range of 0-0.25 GeV. A cut at 0.24 GeV 
was chosen to clean up the E*"*" candidates. This cut pre- 
served ^ 80% of the radiative and 7r° signals, while sub- 
stantially reducing the background under the E*+. Fig. 
[7] shows an example of the Monte Carlo missing energy 
distribution for the jp -> if°E*+ K'^T,+j reaction 
with the dotted line indicating the 0.24 GeV cut. Fig. 
[5] shows the results on the missing mass off the K''^. A 
cut on the invariant mass of the n-7r^ combination along 
with the missing energy cut cleans up the excited-state 
hyperon spectrum, making the E*"^ quite prominent. Fi- 
nally a ±0.03 GeV cut was applied to the missing mass 
off the K'^ around the known mass of the E*+. 

Fig. [5] shows the missing mass squared of all detected 
particles after all of the mentioned cuts. A clear 7r° peak 
is present with some smaller but unknown amount of 
radiative signal at zero missing mass. Fig. 1101 shows the 
the missing mass off the n-7r^ combination. The missing 
mass off the E+ will be used in the background analysis. 



Simulations 

A Monte Carlo simulation of the CLAS detector was 
performed using GEANT [2^, set up for the glla run 
conditions. The experimental photon energy distribu- 
tion for an incident electron beam of 4.0186 GeV was 
used to determine the energies of the incident photons in 
the simulation. Events were generated for the radiative 
channel (E"'"(1385) E+7), the normalization reaction 
(E"'"(1385) E+tt''), and several background reactions. 

A phase space event generator was used with a variable 
t-dependence such that a channel with a K'^ was gener- 
ated uniformly in the center-of-mass frame in (p with a 
t-dependent distribution in 6cm according to P{t) oc e^* 
with 6=2.0 GeV^^. Gaussian distributions in x and y 
with a ~ 0.5 cm were used to approximate the beam 
width at the target. Events were generated uniformly 
along the length of the target. These generated events 
were fed into a simulation of the CLAS detector. 

For each contributing channel the differential cross sec- 
tion was found using data to weight the strength and an- 
gular distribution in the Monte Carlo generator. A very 
careful empirical smearing procedure was used to match 
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FIG. 7: The Monte Carlo missing energy distribution for the 
7P —>■ K^'E*^ — i- K^'E^-f reaction with the dotted line indi- 
cating the 0.24 GeV cut. 
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FIG. 5: The invariant mass of the tt^-tt" (upper left), the 
missing mass off the tt^-tt" (upper right), the u-ttJ invariant 
mass (lower left), and the missing mass squared of all the 
detected particles (lower right). All distributions are before 
any kinematic constraints. 
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FIG. 8: The missing mass off the nf-n with three progres- 
sive cuts applied to isolate the E*^ events. First |A'/(7r^n) — 



< 0.01 GeV, next shown are the events left over after 



the Ex < 0.24 GeV cut, and finally the dotted lines show the 
±0.03 GeV cut around the mass of the E'^. 
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FIG. 6: The invariant mass of the u-ttj with the dotted lines 



showing the \M{tt2 



Mh+I < 0.01 GeV cut (upper left). 



the missing mass off the K'^ (upper right), the total missing 
energy with the dotted line showing the cut used (lower left), 
and the missing mass squared of all the detected particles 
(lower right) after the ±0.01 GeV cut on the peak. 



the Monte Carlo and data resolutions. This procedure is 
discussed in Refs. 



and |27| . Ultimately the missing 



mass squared from Monte Carlo, Mj(i4r°I]+), gave very 
good agreement with the shape of the experimental data 
as shown in Fig. [S] 

This analysis relies on an understanding of the con- 
tributing leakage of background channels into the tt" and 
radiative signal peaks. For example, 7r° leakage from a 
background channel such as 7p — > wA"*" — '!r~ tt^ mr~^ 
will lead to over-counting of — > E+tt^ events. The 
Monte Carlo of various possible contributions was used to 
study the possible background leakage at various stages 
of the analysis. The acceptances of each possible back- 
ground were used to study the possible effect on the final 
reported ratio. 
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FIG. 9: The missing mass squared of all detected parti- 
cles after all analysis cuts. A Gaussian fit gives a mean of 
0.018±0.0002 GeV. 




M,(i:*) (GeV) 
FIG. 10: The missing mass of the E"*". 

At this stage of the analysis the most likely background 
reaction is 7p — > a;A+, followed by A+ mr^ decay. 
The w decays primarily to tt+tt^vt'^ followed by — > 
27. The full reaction jp cl>A+ then has the same 
final state as — s- S+7r° and must be carefully 
considered. This is also true for reactions like jp uN* . 
The A^(1440) has a relatively large decay width at 250- 
450 MeV , which implies possible leakage into any cut 
around the S*"'". This is the reason for the extra steps to 
isolate the E*+ as shown in Fig. [H] Due to the constraints 
on the 7r+ and n, along with the series of cuts shown in 
Fig. m the leakage from the N* decay was negligible. 
However, because the A"*" is very close in mass to the 
S+ and the tt^tt" from the ui decay has a similar phase 
space to the K'^, there was some -fp (u;A+ leakage that 
needed to be accounted for. 

The contribution of the A'*" background was also stud- 
ied. The constraints on the reconstructed K'^, combined 
with the missing mass constraint off the to be the 
mass of the should minimize any K*^ contribu- 



tion. However, because the reaction jp — )■ if*°S"'" has 
the same possible final states that are being analyzed, it 
was carefully considered. The Monte Carlo investigation 
indicated that there were contributions that needed to 
be accounted for. 

Also investigated with Monte Carlo was the reaction 
7P m*N*, where m* is any meson that can decay to 
p7r°, and the N* A^tt"*" provides the detected pion. 
Similarly, •jp — pN*, where the A^* decay to n7r+7r" was 
a possible contaminant. In addition to the kinematic con- 
straints previously mentioned, these backgrounds cannot 
contribute for low W (W < 1.6 GeV). For testing pur- 
poses of these types of reactions, the channel pA^(1520) 
was considered. The p(770) has a width of F = 150.3 
MeV and decays almost 100% to tttt, so it was possible 
to leak under the AT^ invariant mass cut. Ultimately all 
contributions of the channel type jp m*N* were found 
to be negligible (zero acceptance). 

Based on the possible final state decay products, the 
reactions jp TyriTr"*", 7p — )■ if'^E'^Tr"'', and jp — > A'^S"'" 
were also considered. These backgrounds also have 
negligible acceptance as determined from high-statistics 
Monte Carlo studies using the same event selection as for 
the data, and hence were dismissed. 

Minimization of the cj and K*" backgrounds 

As indicated in the previous section, the 7p — > wA+ 
and 7p — ?> Ar*°S+ channels are the most likely back- 
ground contributors. The branching ratio of a; 
7r+7r-7r° is 89.2±0.7% [H], implying a high probability of 
overlap with the normalization channel 7p Ar'^S*^ — > 
7r"'"7r~7r"'"n(7r°). The jp — )■ Ar*°E+ channel was a con- 
cern for the same reason. To get an indication of how 
much these channels were present in the data, the miss- 
ing mass off the 7r^-n combination was used. For the 
7P — ?> a;A+ (7p Ar*''S+) channel the missing mass 
off TT^-n should show a u) (K*^) peak. The missing mass 
spectrum off the tt^-ti combination from the Monte Carlo 
of the 7p — !■ Ar''S*+ — > 7r"'"7r~7r"'"n(7r'^) channel was com- 
pared to the same distribution from data. To isolate the 
tt" channel in the data, a kinematic fit to the missing 
tt'^ with a 10% confidence level cut was applied to leave 
only the final state 7r+7r^7r+n7r" in the data. A direct 
comparison between the data and Monte Carlo of the 
missing mass spectrum off the 7r^-n combination then 
deviated where background was present (normalizing the 
Monte Carlo to the data). It is clear from the comparison 
shown in Fig. [TT] that there is a non-negligible amount 
of oj events. The number of events with a K*'^ present is 
too small to be visible. 

To remove the majority of the uj events a kinematic fit 
was performed with a missing ir'^ hypothesis, while con- 
straining the , tt" and (tt^) to have the ui mass, result- 
ing in a 2-C fit. High confidence level candidates were 
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FIG. 11: The Monte Carlo and data distributions for the 
missing mass off the Trjn combination with the data showing 
a u) peak. The Monte Carlo shown is for the reaction 7p 
K^Tj*^ — >• 7r"*"7r~7r^n7r'', indicating the expected missing mass 
off the TT^n combination if no background was present. The 
dashed lines indicate the known masses of the ui and K**^ . 
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FIG. 12: The Monte Carlo (lines) and data (points) distribu- 
tions for the missing mass off the vr^n combination (as shown 
in Fig. [11]) after the P(x^) < 1% cut showing that the distri- 
butions now match in the uj mass region. 



then rejected as part of the identifiable lu background. 
Various confidence level cuts were tested until the data 
matched the Monte Carlo in the mass range of the a; 
(within the statistical error bars of the data) . Ultimately 
a confidence level cut of P{x^) < 1% was used, resulting 
in the comparison seen in Fig. 1121 The same cut was used 
to reduce the possible 7p ujA'^ 7r+7r~7r+n(7) back- 
ground by imposing the constraint on the 7rj*",7r~ and 
(7) to be uj. In this case it was not possible to use the 
Monte Carlo and data to check in the same way, and so 
^ -P(x^) < cut was used. The same cut was also used 
under a K*'^ hypothesis to reduce the acceptance of the 
7P — > iir*°I]+ — > 7r+7r^7r+ri(7r°) channel. Similarly for 
the "fp — > 7r+7r^7r+n(7) with the hypothesis 
of the 7rj^,7r" and (7) to be K*°. 



Even with the above cuts in place, a small amount of 
the oj and if*" background still slipped through. An es- 
timate of this leakage was found and then subtracted out 
of the final result, as discussed in the following sections. 



Cross sections 

To tune the E*+ Monte Carlo, the differential cross 
sections for the reactions 7p —> K^Y.*^ 7r+7r^ri7r+7r°, 
i<:*°E+ 7r+7r^7r''n7r+, and 7p a;A+ 



IP 



TT ' vr TT mr 



r+ were obtained. The shapes of the differ- 
ential cross sections were then used to adjust the event 
generators. In each case a l/E^ photon energy distribu- 
tion was used in the generator. 

A normalization procedure shown to accurately repro- 
duce a number of well-measured channels was used for 



each cross section [23|. The following is a discussion of 
the procedure used to extract the 7p K^T,*^ cross 
section. A similar procedure was followed for the two 
background channels. 

With all of the aforementioned constraints the jp — )■ 
^Oj^*+ j{0-^+^o reaction was easily isolated with a 
kinematic fit to the missing tt*^. A 10% confidence level 
cut was applied to ensure channel purity. The yield was 
determined by the ratio of the raw S+tt" events to the 
number of incident photons in each E-^ bin, so as to nor- 
malize with the bremsstrahlung spectrum. Corrections 
were made for each bin with the newly obtained accep- 
tances. The angular dependence in the generator was 
initially flat with a zero i-slope dependence. After the 
differential cross section was obtained, the distributions 
were used to adjust the S*+ generator. Each correspond- 
ing angle and energy bin was filled according to the dis- 
tributions seen in the data results. Each angle bin was 
divided into E-y bins and represented accordingly in the 
new event weighting scheme of the generator. The ad- 
justed generator was then used to produce new Monte 
Carlo and obtain more accurate acceptances. This pro- 
cess was then iterated until no change was seen to the 
differential cross sections within the statistical uncertain- 
ties. After these modifications were made, the resulting 
Monte Carlo was compared with the data, using the mo- 
mentum distributions for the vr^, tt^, and neutron tracks, 
as well as the K° lab frame angle distribution, and found 
to closely match the data within the statistical uncer- 
tainties, see Fig. [T3l 

The same corrections were applied to the jp — > 
K°'E,*+ -)• K°T,+-f cross sections. The nature of the 
corrections to the Monte Carlo were specific to the jp 
iir°I]*+ cross section and so the corrections could be ap- 
plied without discrimination between the electromagnetic 
and strong decays of the baryon. 

To calculate the acceptance of the signal and back- 
ground reactions, an extraction method used to resolve 
the radiative and tt" channels was required and will be 
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FIG. 13: The momentum distributions for the vr'*', 7r~, and 
neutron are shown in the first 3 panels, with the cosine of the 
kaon lab frame angular distribution in the last panel, for the 
reaction jp — > A'^E*^ -> tt' tv'^ mr^ tt'^ . The data and Monte 
Carlo are shown as the histogram and points, respectively, 
which closely overlap. 



discussed next. 



Fitting technique 

The two-step kinematic fitting procedure developed in 
Ref. [l^ was employed to resolve the radiative and strong 
decay signals. Because of such similar topologies and the 
small relative size of the radiative signal, the kinematic 
fitting procedure could not be expected to cleanly sepa- 
rate the E"^7 events from the overwhelming S+tt'^ events 
in a single fit. Thus we employed a two-step kinematic 
fitting procedure, making first a kinematic fit to a miss- 
ing 7r° hypothesis, then checking the quality of the fit of 
the low confidence level candidates in a second kinematic 
fit to the actual radiative hypothesis. 

In order to check the quality of the kinematic fit to 
a particular hypothesis, we studied the distribution 
from the fitting results. In this procedure all detected 
particles were kinematically fit to the appropriate miss- 
ing mass hypothesis. An additional constraint was in- 
troduced into the kinematic fitting enabling analysis of 
the more well-behaved 2-C distribution, as opposed 
to the 1-C distribution, to test the quality of the candi- 
dates with the hypothesis used. The constraint required 
that the neutron and 7r+ track have an invariant mass of 
the E"*" in the hypothesis. The detected particle tracks 
were kinematically fit as the final stage of analysis and fil- 



tered with the confidence level cut. In this fit there were 
three unknowns (px) and five constraint equations, four 
from conservation of momentum and then the additional 
invariant mass condition. This makes a 2-C kinematic 
fit. In the attempt to separate the various contributions 
of the E*"'' radiative decay and the decay to S+7r°, the 
events were fit using different hypotheses for the topol- 
ogy: 

7P — > 7r+7r^7r+n(7r'') 2-C 
7P — > 7r+7r^7r+n(7) 2-C. 



The constraint equations were 



beam ~t~ ^'-^p 
Pbeam 



Et^+ - E^- — — En — Ex 

-P7T+ -PtT- -PtT+ ~Px 



0. 



(5) 



Px and Ej. represent the missing momentum and energy 
of the undetected 7r° or 7. 

Then a fit function was made for a 2-C distribution 
following from Ref. ^] as 



fix') 



2 



(6) 



This fit function has a fiat background term, P2. Pi 
was used to measure how close the distribution in the 
histogram was to the ideal theoretical distribution for 
two degrees of freedom. 

Because there were two kinematic fits for both the 7r° 
and radiative channels, some new notation is introduced. 
The first confidence level cut used to filter out the larger 
tt" signal from the radiative signal by using a kinematic 
fit to I]"'"(7r°) and taking only the low confidence level 
candidates is denoted as P^(x^)- The final kinematic fit 
used to isolate the radiative signal, using a S+(7) hy- 
pothesis has a confidence level cut denoted as P^(x^), 
taking only the high confidence level candidates. Opti- 
mization studies have been previously done to constrain 
the choice of P^(x^) and Pi;(x^) [26|. 



Ratio calculation 

The TT*^ leakage into the 7 channel was the dominant 
correction to the radiative branching ratio. To prop- 
erly calculate the ratio, the leakage into the tt'^ region 
from the 7 channel was also used. Taking just these two 
channels into consideration, the number of true counts is 
represented as N{Tj^) for the S*"*" — )■ S+7 channel and 
N{T,tt) for the S]*+ — > S+7r° channel. The acceptance 
under the S*"*" E+7 hypothesis is written as A^{X), 
with the subscript showing the hypothesis type and the 
actual channel of Monte Carlo input that was used to 
obtain the acceptance value indicated in the parenthe- 
ses. For the calculated acceptance of the E*+ — >■ E+7 
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channel under the S*'^ — > S+7 hypothesis, the accep- 
tance is A^CEj), and for the !]*+ — > I]+7r° hypothesis it 
is A^{'Ej). It is now possible to express the measured 
yields for each channel and 

= A^(E7)7V(E7) + A^{J:tt)N{T.tt) (7) 

= A,(S7r)Ar(E7r) + Ar,iEj)N{i:-f). (8) 

The desired branching ratio of the radiative channel 
to the 7r° channel using the true counts is then R = 
N {H'j) / N (Y^tt) . Solving for R to get the branching ratio 
expressed in terms of measured values and acceptances 
gives 

^ _ n^A^{ETr) - n^A^jEn) 
n^A-,(T,-f) - n^A^CEj) ' 

Equation [9] is based on the assumption that there are no 
further background contributions. The formula for the 
branching ratio to take into account background from 
the cjA"*" channel, as an example, can be expressed as 

^ An-yAT^(T,Tr) — AriTrA^CET:) 

(10) 

where 

- N^{uj n+n-j) (11) 

and 

Arij = n,y — N^{uj — > 7r~'"7r~7) 

- N^{u} 7r+7r"7r°). (12) 

The n-y (tItt) terms come directly from the yield of the 
kinematic fits and represent the measured number of pho- 
ton (pion) candidates. In the notation used, lowercase n 
represents the measured counts, while uppercase N rep- 
resents the acceptance corrected or derived quantities. 
The A^^.TT terms are corrections needed for the leakage 
from the wA+ channel (an arbitrary number of back- 
ground types N-y^T^{X) can be accounted for in this man- 
ner). The notation utilized is such that the pion (photon) 
contributions are denoted N^^ (-^7), so that Nj{X) de- 
notes the relative leakage of the {X) channel under the 
E+7 hypothesis and N^^iX) denotes the relative leakage 
of the X channel under the S+tt^ hypothesis. 

The final acceptance for each channel was determined 
after the final set of confidence level cuts was taken. Af- 
ter the background acceptances were minimized, an esti- 
mate of the background contributions was found for each 
relevant case. 



TABLE I: The cuts used to extract the final radiative and tv' 
counts. (See text for details.) 



Cut Used (Appfied) 

(1) |Af(7r+7r-) ~ Mko \ < 0.01 GeV (both) 

(2) |M(7r+n) - \ < 0.01 GeV (both) 

(3) < 0.24 GeV (both) 

(4) |M,(7r+7r-) - M^,+ \ < 0.03 GeV {both) 

(5) P^"(x') < 0.01% {7v°) 

(6) P^^ix') > 10% (^") 

(7) P;(x') < 0.01% (7) 

(8) R%x') > 10% (7) 



SIGNAL EXTRACTION 

Each Monte Carlo channel was run through the analy- 
sis with the same cuts as used for the data. These cuts for 
the extraction of the radiative and 7r° signals are listed 
m Table H The cuts are listed in the order implemented. 
The first cut, (1), was on the mass and restricted the 
TTj^, TT^ sample. The second cut, (2), on the S+ mass was 
implemented to clean up the tt^ , n sample prior to the 
more restrictive cuts, (5)-(8). The missing energy restric- 
tion used to reduce background is number (3). Cut (4) 
restricted the missing mass off the K'^ to be in the range 
of the Cuts (5) and (6) list the final confidence level 
cuts used from the kinematic fit to isolate the missing tt'^. 
Cuts (7) and (8) were used to isolate the radiative decay. 
The second column lists whether the cut was applied to 
just one channel or both. 

The two-step kinematic fitting procedure was used to 
isolate the radiative signal from the 7r° channel. In this 
procedure two separate kinematic fits were performed, 
one with zero missing mass for the 7 and the other with 
the missing mass of the 7r°. The fit function in Eq. [6] was 
used to fit the distributions to determine the resulting 
quality of candidates present in the fit. The parameter 
Pi was used to measure how close the distribution in the 
histogram was to the ideal theoretical distribution for 
two degrees of freedom. The pure radiative decay Monte 
Carlo was used to determine the value of the expected 
Pi parameter. Fig. [T3] (A) shows the distribution 
for a kinematic fit of the Monte Carlo channel jp 
^oj]*+ _^ ^+^-^+^^0 under the radiative hypothesis, 
displaying a highly distorted distribution. Fig. [14] (B) 
shows the x^ distribution from the same kinematic fit 
under a radiative hypothesis of the Monte Carlo channel 
7P —5- Ar°S]*+ 7r+7r^7r+n7, a fit using Eq. [H] after all 
the kinematic cuts up to (4) in TableU The parameter Pi 
was used as the ideal value to expect. After obtaining the 
expected Pi, a kinematic fit to the data was performed 
using both hypotheses. 

The first confidence level cut P^{x^) was used to filter 
out the larger tt^ signal from the radiative signal by using 
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a kinematic fit to a I]+(7r°) hypothesis and taking only 
the low confidence level candidates. This confidence level 
cut was checked and made more restrictive until Pi from 
the data matched the expected value from Monte Carlo 
(within the projected error bars). Fig. [T3] (C) shows 
the distribution and fit before any P^(x^) cut was 
applied and Fig. [14] (D) shows the distribution after a 
P^hd) < 0.01% cut was applied. The final confidence 
level cut i-'^(x^) from the kinematic fit to a S+(7) hy- 
pothesis was used on the remaining candidates. Only the 
high confidence level candidates were preserved. 

Note that the yields for the S+tt'^ decay will be re- 
duced for a lower value of P^(x^), which is desirable for 
extracting the radiative decay signal. On the other hand, 
this cut cannot be made arbitrarily small, since it reduces 
the statistics {i.e., increases the statistical uncertainty). 
Similarly, the E"*'7 signal will be purified by a higher 
cut on P^(x^), but again the higher the cut, the lower 
the statistics. The Monte Carlo was used to examine 
the acceptance of these cuts for various branching ra- 
tios (E+7/S+7r°), which is discussed in the next section. 
Ultimately, the branching ratio extracted from the data 
should not depend on the cut points chosen (assuming the 
Monte Carlo gives accurate cut acceptances). The Monte 
Carlo was then used to optimize the trade-off between 
statistical uncertainty and systematic uncertainty (due to 
the choice of confidence level cuts based on a more quan- 
titative analysis of Pi). The cut value of P"(x^) < 0.01% 
showed consistent optimization with P^{x^) > 10% for 
this topology and range of statistics as listed in Table 
|I1 Details of the optimization method of the confidence 
level cuts using the Monte Carlo are described in Ref. 
[2^ . The confidence level distribution under the radia- 
tive hypothesis before and after the P^(x^) < 0.01% cut 
is shown in Fig. [TSjA). Likewise, the confidence level 
distribution for a fit to the data under the missing tt*^ 
hypothesis before and after the Pi'(x^) < 0.01% cut is 



shown in Fig. [TSiB). 

The same cuts determined to effectively isolate the 
radiative signal, P^(x^) < 0.01% with P^(x^) > 10%, 



were used to isolate the 



channel used for normal- 



ization of the ratio, such that P"(x^) < 0.01% is used 
with P^(x^) > 10%. This can be denoted as P^{x'^) = 
Pjix^) < 0.01% and P^ix^) = P^ix^) > 10%. The final 
missing mass squared distribution after all cuts is shown 
in Fig. [TBI before the two sets of confidence level cuts and 
after. 

The acceptances were found for each contributing 
channel and are listed in Table [iTj Each of the chan- 
nels in Table [U was generated with enough statistics so 
that the statistical uncertainty would not contribute in 
the final ratio calculation. The value of the acceptance 
for the 7r° (7) hypothesis is listed under the (Aj) 
column. The uncertainty is statistical only. 

The acceptance values indicate that contributions from 
the and uj channels will be subtracted out directly. 
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FIG. 14: Plot (A) shows the distribution for the Monte 



Carlo channel 7p —>■ E* 



■ TT TT TT nyr 



under the radia- 



tive hypothesis, displaying a highly distorted x distribution. 
Plot (B) shows the fit results for "/p — >■ K°Tj*'^ -k^ -k^ wy , 
displaying a reasonable distribution with the ideal Pi . The 
radiative hypothesis kinematic fit distribution and distri- 
bution fit to data before any P"(x^) cut is applied is shown 
in (C) and the same distribution after a P^ix^) < 0.01% cut 
is applied is shown in (D). 




0.2 0.4 0.6 0.8 1 

P(X') 



0.2 0.4 0.6 0.8 1 



FIG. 15: (A) The confidence level distribution for a fit to 
the data under the radiative hypothesis before the -P"(x^) < 
0.01% cut (line) and after (dotted line). (B) The confidence 
level distribution for a fit to the data under the missing tt" 
hypothesis before the P"(x^) < 0.01% cut (line) and after 
(dotted line). 



All other background channels not listed were zero. As 
mentioned, all background contributions to the ratio are 
relatively small, but care is taken to accurately consider 
each contribution. The levels of these contributions de- 
pend the placement of the confidence level cuts previously 
mentioned. 

To obtain an estimate of the amount of leakage into the 
S"'"7 and S"'"7r'^ signals, some cuts were removed to obtain 
a fit on the channels of interest. Only the |M(7r||''7r^) — 
Mko \ < 0.01 GeV and the |M(7r+n) - Ms+ 1 < 0.01 GeV 
cuts from Table[T]were used with an additional cut on the 
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FIG. 16: The missing mass ofT of all detected particles. The 
plot shows the final radiative candidates at zero missing mass 
after the P"(x^) < 0.01% and P^{x^) > 10% cuts. Also 
shown are the final vr" candidates after the -P"(x^) < 0.01% 
and P'(x^) > 10% cuts. 




M,(7t*n) (GeV) 

FIG. 17: The missing mass off the vrj-n combination with 
the \M{ntn-) - M^o \ < 0.01 GeV, |Af(7rf n) - Mh+ | < 0.01 
GeV, and \Ml - M^o| < 0.0175 GeV^ cuts. The resuhing 
K*^ and uj peaks are fit with Breit-Wigner line shapes, while 
the background is fit with a polynomial function. 



TABLE II: Acceptances (in units of 10~^) for the channels 
that survive all cuts. All of the cuts used to obtain the ac- 
ceptance values are listed in Table HI The uncertainties are 
statistical only. The two columns contain the acceptance for 
each hypothesis At^, A-y. 

Reaction ^4^ A^, 

JC"E*+ ^ A'"E+7 0.0644±0.0040 0.6244±0.0125 

K°T.*+ -> K°Y.+-n° 0.6502±0.0128 0.0591±0.0038 

is:*°E+ K°T.+"i 0.0018±0.0005 0.0186±0.0021 

iC*°E+ ^ K°Y.+-n° 0.0231±0.0023 0.0050±0.0011 

tjA+ -> 7r+7r"7r°n7r+ 0.0003±0.0001 O.OOOOiO.OOOO 

a;A+ 7r+7r~7W7r+ 0.0000±0.0000 0.0002±0.0000 



missing mass squared of all detected particles around the 
7r° mass of jM^-M^ol < 0.0175 GeV^. The missing mass 
off the TT^fT. combination was then checked. The resulting 
K*^ and w peaks were fit with a relativistic Breit-Wigner 
while the background was fit with a polynomial function, 
as shown in Fig. [T71 

The total number of K*^ events present in the data set 
using the less restricted set of cuts just described can be 
expressed as 

N{K*^) = 

R{K*^ TT+'K-TT°)A^'\K*° tt+tt-tt") ~ 

3.639 X 10*^ ± 6.6 X 10"*, (13) 

where n{K*° tt+tt-tt") = 3019 ± 55 is the esti- 
mated number of 7p — > K*^Yi'^ Tr+Tr^vr^Tr'^n events 
found through the integrated fit to the K*^ in Fig. [T71 
R{K*^ 7r+7r-7r°) = 1/3 is the probability for the de- 
cay and A^ (A'*" — 7r+7r~7r'') is the acceptance of the 
K*^ — > TT+TT^TT*' channel found by observing how many 



thrown events survive the three cuts used to obtain the 
K*^ sample. The value N{K*°) can then be used to ob- 
tain an estimate of the K*'^ contribution from any set 
of cuts, given an accurate acceptance for the new cuts. 
The number of K*'^ events that would be present in the 
analysis outlined in Table U under the tt^ hypothesis can 
be expressed as 

N^{K*" tt+tt-tt") = 

R{K*° ^ tt+t:-tt°)A^"'{K*^ ^ tt+tt-ttO) ' 

where {K*'^ — > tt+tt^tt'^) is the acceptance for the 
7P A:*°S+ 7r+7r^7r"7r+n channel under the S'''(7r°) 
hypothesis. Likewise for the ^^(7) hypothesis inter- 
changing A^{K*^ 7r+7r-7r°) with A^{K*° tt+tt-ttO) 
to obtain N^{K*° 71+71-71°). The radiative decay 
can also have a contribution under the I]+(7r°) hypothe- 
sis, 

R{K*° 7T+7T-j)A^{K*° 7T+7T-7T°)N{K*"), (15) 

or under the ^+(7) hypothesis, 

N.f{K*" 7r+7r"7) = 
R{K*° ^ 7r+7r--f)A^{K*" ^ 7r+7:- 7r")N {!<*"), (16) 

where R{K*° 7r+7T--/) = 2.39 ± 0.21 x lO'^ |2^]. 

In the case of the lo contributions, no distinction is 
made between the jp — )■ ujA+ and jp — )■ ujn7:+ chan- 
nels. The Monte Carlo used in the background esti- 
mate for the uj is the jp — wA"*" channel only because 
there is a slightly larger acceptance for this channel. 
By using the channel of greatest acceptance an over- 
estimate is expected. The total number of events from 
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wA"*" — ;> TT+TT 7r°n7r+ present is estimated as 

N{uj) = 

n{uj vr+TT^vr'') 

5.296 X 10^ ±4.7 X 10^ (17) 

Here n(a; — >■ 7r+7r^7r°) = 11120 ± 106 is the estimate 
from integrating the w fit in Fig. [T71 R{uj 7r+7r~7r°) = 
89.2 ± 0.7% is the branching ratio of the w decay to 
7r+7r"7r° Q and -> 7r+7r"7r°) is the probability 

that this decay channel will be observed after the three 
cuts used to obtain the fit to the w peak. An estimate 
of the number of counts under the S+(7r°) hypothesis 
coming from the uj is obtained using Eq. [17] as 

A^{uj 7r+7r-7r")i?(w 7r+7r"7r°)A^(w), (18) 

where A^(a; ^> tt+tt^tt^) is the acceptance for the 7p — > 
a;A+ — !> 7r+7r^7r'^7r+n channel under the E"'"(7r'^) hypoth- 
esis. 

It is possible to express all other associated lo correc- 
tions to be obtained given the value of iV^, along with 
the acceptance terms for that particular channel. The 
corrections for the 7 and tt channels, respectively, are 
written as 

iV^,^(a; ^ tt+tt-ttV) = 
A^^^{uj 'k+ti-'k"'k^)R{uj 7r+7r"7r°7r°)A^(w), 

N^i^T^{uj — 7r^7r^7) = 
A^^T,{uj 7r+7r^7)i?(a; ^ 71^11^ ^)N{uj), 

N^^T,{uj tt+tt^tt") = 
A^^^{u} TT+Tr-TT°)R{uj n+Tr-7r")N{u}). (19) 

where R is used for the corresponding branching ratio or 
upper limit in each case, for example, R{uj — > 7r"'"7r~7) 
is the branching ratio for the radiative decay of the lu 
with a value less than 3.6 x 10"'^. The value of R{uj 
7r+7r-7r°7r") is listed at less than 2 x 10"'* ^24]. All results 
from background contributions are tabulated in Table lTTll 
All non-listed background is considered negligible. 

Final yields 

To calculate the ratio of the EM decay to the strong 
decay, Eq. |9] is employed. All terms that take into ac- 
count any channel other than the and radiative signals 
are for the time being ignored. The acceptance values are 
taken from Table HT] The raw values obtained out of the 
final kinematic fit are = 148 and = 682, as seen 
in Fig. 1161 with statistical uncertainties taken as the 
square root of n in each case. After accounting for the 
backgrounds listed in Table IIIH the corrected counts are 
An^ = 135.81 ± 11.99 and An^ = 642.11 ± 26.38. 



The ratio of the K°T,*+ K°T,+-f channel to the 
K°J:*+ K"j:+Tr° channel is then, 

nS*+^s+7 _ An^A^(i:7r) ~ An^A^{T,Tr) 

= 11.95 ±2.21%. (20) 

The raw counts for the radiative and tt'^ extraction 
were obtained using ~ < 0.01% with the final 
confidence level cut P}^ ~ P^ > 10% as mentioned in the 
cuts from Table HI 

Only the statistical uncertainty is quoted. To deter- 
mine how reliable the ratio is a set of systematic studies 
is required along with a study of the variation in the ratio 
based on the choice of confidence level cuts. This vari- 
ation and all other systematic studies are considered in 
the next section. 

SYSTEMATIC STUDIES 

The value of each of the nominal cuts was varied to 
study the effect on the final background corrected ratio. 
For each variation the new acceptance terms in Equation 
iniwere recalculated with the corresponding Monte Carlo. 
Each major systematic uncertainty contribution is num- 
bered as it is discussed and listed in Table IVIIIl which 
contains a summary of all systematic uncertainties. 

Several A/3 cut variations were checked starting with 
A/3| < 0.02 for all charged particles, leading to a ratio 
of 11.98±2.22%. There was also a check at |A^| < 0.1 
that gave a ratio of 11.74±2.17%. The |A/3| cut selected 
uses a ±1 ns timing cut, while keeping |A/3| < 0.035 for 
all pions. This variation is presented in Table IVIIIl as 
number (1). 

To estimate the systematic effects from the Monte 
Carlo, such as the uncertainty in correctly simulating the 
data, a comparison was made with the cross section of 
7P — !■ K^Y.*^ from Monte Carlo and data. The ratio 
was obtained using the acceptance corections based on a 
Monte Carlo with a zero i-slope in the generator to get 
R = 11. 70 ±2. 21%. This value deviated from the ratio ob- 
tained by ~ 2%. The uncertainty from Monte Carlo was 
then estimated to be ~ 2% in either direction. This find- 
ing leads to a high value in the ratio of i? = 12.19±2.21% 
and a lower value in the ratio of i? = 11.71±2.21%, which 
is listed as number (2) in Table Ivnll 

The missing energy cut removes a large amount of 
background that slips in under the E*"*". The final cut 
of Ex < 0.24 GeV was chosen based on the maximiza- 
tion of signal counts (decreased statistical uncertainty). 
The systematic contribution for the cut was stud- 
ied by varying the cut in a reasonable range about the 
nominal value. This ratio was reasonably stable up un- 
til Ex < 0.28 GeV, at which point the 7r° and 7 sig- 
nals become overwhelmed with background. Cuts in Ex 
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TABLE III: The contributions for all of the background channels taken into consideration. The uncertainties are statistical 
only. The two columns contain the contributions estimated for each hypothesis A^.^, A'^^. 

Reaction A^,r N-y 

K"'T.+ i^"E+7r" 28.02±2.84 6.02±1.34 

K*°T.+ A'°E+7 0.0157±0.0046 0.162±0.023 

a;A+ ^ 7r+7r-7r%7r+ 11.81 ± 2.37 0.0120 ± 0.0024 

ujA+ Tv+n-'ymv+ 0.0400 ± 0.0057 2.1 x 10"^ ± 3.0 x 10"'' 

LjA+ TV + TT-TV°TT°mV+ 1.1 X 10"* ± 2.0 X 10"^ 0.0 



too low in energy tend to distort the ratio. Based on 
this study a high (11.95±2.21%) and low (11.71±2.20%) 
value were assigned to the associated systematic uncer- 
tainties. These contributions are listed in Table IVllIl as 
number (3). 

The background counts that contribute to the ratio 
assume a branching ratio for u — 7r"'"7r~7 to be at 
the top of the upper limit at 3.6 x 10~^, as well as for 
" at 2 X 10~^ [2J|. The variation in uncer- 



tainty of the branching ratio of either of these channels 
is not of large enough order to make a notable effect on 
the ratio. The total contribution from background in the 
7r° hypothesis was 39.89 ± 3.70 and for the 7 hypothesis 
was 6.19 ± 1.34. To check for the largest deviations under 
these uncertainties, a contribution of 43.59 counts from 
the 7r° hypothesis was used with 4.85 counts for the 7 
hypothesis, leading to a ratio of 12.37 ± 2.24%. The op- 
posite extreme was also used to obtain 36.19 counts from 
the tt'^ hypothesis and 7.5 counts for the 7 hypothesis, 
leading to a ratio 11.54 ± 2.19%. This finding is listed in 
Table Em] as number (4). 

The kinematic fits used to control the leakage of the 
K*'^ and oj have an associated confidence level cut that 
was also tested. The cut in each case was selected based 
on the reduction of the specific background channel, 
while maximizing the signal counts from the 7 and tt'^. 
A check was done over a large range of confidence level 
cuts for each case to test the variation in the ratio. New 
background contributions were found for each cut along 
with new acceptance terms. The ratio was then recalcu- 
lated and tabulated for the K*'^ (w) removal in Table ITVl 
(Table \V\ . These contributions are listed as (5) and (6), 
respectively, in Table IVlTll 

The optimum set of confidence level cuts used to ex- 
tract the final yields has a range of validity in which a 
study of the ratio variation is appropriate. The range 
of validity is found by using the fractional deviation in 
the ratio 6R, and requiring it to be less than or equal to 
the fractional uncertainty due to statistics. SR is defined 
as the difference in the generated ratio and recovered ra- 
tio in the Monte Carlo study. The set of optimum cuts 
occurs at different values for a given mixture of Monte 
Carlo S+TT^ and S"'"7 events. Using the glla data, the 
Monte Carlo was tuned to have approximately the same 





R 


0.100 


12.04±2.20% 


0.050 


11.79±2.20% 


0.010 


11.95±2.21% 


0.005 


11.92±2.22% 


0.001 


11.90±2.21% 



TABLE IV: The values used in the systematic variation for 
the calculated ratio with respect to the confidence level cut 
from the fit to jp — >■ iv'^nK*'^ — >■ 7r^n7r"'"7r~ (vr'') and 7p — >■ 
Tv'^nK*'^ — > 7r+n7r+7r" (7). The confidence level cuts were set 
equal to each other such that P-yix^) = fV(x^)- 



R 



0.100 


11.99±2.20% 


0.050 


12.01±2.21% 


0.010 


11.95±2.21% 


0.005 


11.80±2.22% 


0.001 


11.77±2.21% 



TABLE V: The values used in the systematic variation for the 
calculated ratio with respect to the confidence level cut from 
the fit to 7p — >■ Tv'^nui Tv'^mT^-K~ {tv") and yp -K^nuj 
7r"'"n7r+7r~ (7). The confidence level cuts were set equal to 
each other such that P-yix?) = Pnix^)- 



ratio and the same statistics as the real data, and the op- 
timum cuts were thereby determined quantitatively. This 
set of cuts used in the systematic studies, along with the 
values of are given in Table IVII for isolation of the 
radiative decay and Table IVIII for the tt" isolation. Note 
that the value of 5R is not the systematic uncertainty in 
R^-ther, the systematic uncertainty comes from 
the data ratio values in the right-most column in each 
table. The systematic uncertainty is based on the vari- 
ation in the extracted ratio of Eq. (fTO)) for a set of cuts 
determined to give the minimal deviation in the Monte 
Carlo recovered ratio for the valid range of statistics. 

The systematic dependence on the choice of the P[x^) 
cuts in both the radiative and 7r° hypotheses comes from 
Table and Table lylll The ratio was recalculated from 
the resulting raw counts in each case with the new ac- 
ceptance terms for the set of cuts obtained from the op- 
timal range. The set of cuts was tested for both the ra- 
diative and 7r° hypotheses separately. In each variation 
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R 



0.050 17 0.140 11.68±2.22% 

0.005 15 0.090 11.88±2.29% 

0.010 10 0.089 11.95±2.21% 

0.075 20 0.092 11.41±2.27% 



TABLE VI: Optimization points for each P° and P!^ for the 
E"'"7 channeL 



Py{%) P^{%) SR 



R 



0.050 17 0.187 11.35±2.23% 

0.005 15 0.109 12.13±2.31% 

0.010 10 0.099 11.95±2.21% 

0.075 20 0.195 11.19±2.30% 



TABLE VIL Optimization points for each P° and for the 
E^TT*^ channeL 



under the 7r° (7) hypothesis, the opposing cut for the 7 
(7r°) hypothesis was not changed. The highest and lowest 
value from each study were used as the contributions to 
the systematic uncertainty hsted as (7) and (8) in Table 

Eull 

The cut on the tti-tt^ invariant mass was used to min- 
imize the events that are not associated with a K^. A 
range of cuts on \M{Tr^n~) — MKo\ was used to study this 
effect. Only the stable region was used in the final deter- 
mination of the range of variation for this cut with a high 
value of 11.95±2.21% and a low value of 11.72±2.21%. 
These are hsted in line (9) of Table EnH 

The cut on the tt^-ti invariant mass regulates the can- 
didates going into the final set of kinematic fits so there 
is some sensitivity to any events that do not come from 
the E+. A range of cuts on |Af(7r^7i) — M^+\ was used 
to study this effect. Only the stable region was used in 
the final determination of the range of variation for this 
cut with a high value of 11.95±2.21% and a low value of 
11.34±2.85%. These are listed in line (10) of TableEnH 



Table IVini shows a summary of the systematic studies 
and the highest and lowest value of the ratio based on 
the variations mentioned for each type of uncertainty. To 
calculate the final systematic uncertainty, the difference 
in the ratio R = 11.95% and the high value of the ratio 
for each case in Table IVIIII was added in quadrature to 
obtain a value for the uncertainty of 0.53% greater than 
the ratio. The lower systematic uncertainty bound was 
based on the difference between the ratio R = 11.95% 
and the low value of the ratio for each case, resulting 
in a value of 1.21% less than the ratio. The final ratio 
reported is 11.95 ± 2.21(siai)t?;2?(s2/s)%. 



OVERALL RESULTS 

The final result for the ratio of the S*+ - 
E*+ —5- S+TT*^ with systematic uncertainties is 



i;+7 to 



Rt 



11.95 ± 2.21{stat)tifi{sys)%. (21) 



To calculate the EM decay partial width from the mea- 
sured branching ratio, the full width of the E*"*" decay is 
used, rj;^„/, = 35.8 ± 0.8 MeV [H, with the branching 
ratio i?(E*+ ^ E+ttO) = 5.85±0.75%. The partial width 
calculation including systematic uncertainties leads to, 



r 



R 



loR{^* 



250±56.9{stat)tl'^iisys) keV. (22) 

To obtain the corresponding U-spin prediction, we first 
look at the prediction for the A+ p-f partial width to 
the E*+ ^ E+7 partial width. 



(A+|p7) = 
(E*+|E+7) = 
leading to a ratio of, 

r(A+ ^ p7) 



2'2 

i|i 
2'2 



1 

- X 



loo 



= 1 
= 1, 



^ [Mil 

S+7) \Ma 



= 2.638. 



The value for the center of mass momentum for the pro- 
ton is qp = 0.259 GeV and for the E+ is q+ = 0.180 GeV, 
23. 

This implies that the U-spin prediction for the partial 
width of the electromagnetic decay using the 660 ± 60 
MeV width of the A+ decay [2J| is. 



r(A+ ^ P7) 

2.638 



250 ± 23 keV. 



(23) 



A similar calculation can be done for to obtain a U-spin 
prediction for the T-^to_^\^ partial decay width. Table 
IIXI shows the previous model predictions along with the 
U-spin prediction and the final results from this analysis 
in each case. The statistical and systematic uncertainties 
are combined in the present partial widths. 

The partial width in both cases is larger than any pre- 
diction listed except for the U-spin prediction. The U- 
spin predictions for the ^ A7 and E*+ — > E+7 par- 
tial widths are well validated by the experimental result. 
For these cases U-spin symmetry is confirmed within the 
experimental uncertainties. It is important to note that 
the U-spin prediction for the E*"*" EM decay partial width 
ignores the effects of the interference of the isovector and 
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TABLE VIII: Sources of systematic variation in the ratio (in %) showing the contributions to the systematic uncertainties from 
changes to the event selection values, along with the low and high value used. 



Source 


Low Value 


Low Contribution High Value Hif 


Contribution 


(1) A/3 


11.74±2.17 


-0.21 


11.98±2.22 


-K0.03 


(2) MC-match 


11.71±2.21 


-0.24 


12.19±2.21 


-hO.24 


(3) 


11.71±2.20 


-0.24 


11.95±2.21 


-HO.OO 


(4) BG uncertainty 11.54±2.19 


-0.41 


12.37±2.24 


4-0.42 


(5) A'*°-CL 


11.79±2.20 


-0.16 


12.04±2.20 


-h0.09 


(6) uj-CL 


11.77±2.21 


-0.18 


12.01±2.21 


-H0.06 


(7) 7-CL 


11.41±2.27 


-0.54 


11.95±2.21 


+0.00 


(8) 7r°-CL 


11.35±2.23 


-0.60 


12.13±2.31 


+0.18 


(9) K'^ cut 


11.72±2.21 


-0.23 


11.95±2.21 


+0.00 


(10) E+ cut 


11.34±2.85 


-0.61 


11.95±2.21 


+0.00 


Total Uncertainty 




-1.21 




+0.53 



TABLE IX: Comparison of theoretical model predictions for 
the radiative decay widths with the experimental result for 
E*" from Ref. [l^ and the present result result for E*"*". 



Model 


E(1385)'' A7 E^ 


-(1385) -> E+7 


NRQM [3, 4, 7| 


273 


104 


RCQM [5] ~ 


267 




XCQM [6] 


265 


105 


MIT Bag [7] 


152 


117 


Soliton [8| 


243 


91 


Skyrme [9, 10\ 


157-209 


47 


Algebraic model [11] 


221.3 


140.7 


U-spin 


423+38 


250+23 


Results 


445+102 


250+70 



isoscalar components of the photon. If the isoscalar com- 
ponent interfered destructively, the resulting prediction 
could indeed be much smaller. 

The results in Ref. reveal that the meson cloud 
effect can contribute significantly (~ 40%) to the overall 
electromagnetic decay width of the A Nj. This puts 
the prediction from the model at about 80% of the ex- 
perimental measurement. As stated previously it has not 
yet been determined from a theoretical standpoint if the 
meson cloud effects contribute and if so to what degree 
for the radiative decay of the E*" and E*"*". This may be 
the reason for such a difference in the predictions seen 
from experiment compared to the models listed in Table 
HXl Because the U-spin prediction for the E* EM decay 
width uses empirical information from the A EM decay, 
contributions from phenomena like the meson cloud ef- 
fect should be inherent. The correspondence between the 
U-spin prediction and the experimental result strongly 
suggests that the other models lack that corresponding 
effect. 

Perhaps this work can prompt more encompassing cal- 
culations that are necessary to probe the structure of 
the baryon resonances and motivate consideration of the 



meson cloud contributions for electromagnetic decay pre- 
dictions. 
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